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SCATTERING from condensates in turbulent jets 


I. INTRODUCTION 

Measurements ^ ^of turbulent intensities in supersonic 
jets, using an ultra-violet, crossed-beam system appeared to 
give a signal which originated in scattering of light from 
condensed water vapor or some other condensate, rather than in 
the absorption of light by 0 2 molecules. Further measurements 
with a He-Ne laser crossed-beam instrument^ 2 -* substantiated 
this observation. It was in part this uncertainty about the origin 
of the signals, that led to the development of an infra-red 
crossed-beam system. The infra-red system has the advantage 
of the ultra-violet system in that it is sensitive to a 
naturally- occuring air specie (C0 2 rather than C> 2 ) but also it 
operates at a more favorable wavelength for the suppression of 
scattering. Whether or not a real gain is realized, however, 
depends upon the relative levels of scattering and absorption 
signals. It is the purpose of the study reported here to 
develop a mathematical model for scattering from condensates 
and to apply the results to assess the contribution of 
scattering from condensed water droplets in the wavelength 
region of the C0 2 (4.3p) vibration band. 
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II. DEVELOPMENT OF THE SCATTERING MODEL 


A. Scattering Processes 

Light scattering from particles results from the 
interaction between the electric field of the light wave and 
the electric dipole moment (either permanent or induced) in 
the particle. The effectiveness of a particle as a scatterer 
is described by its scattering cross-section K, which is the 
ratio of the effective cross-sectional area of the particle 
in reducing the incident light intensity, to its actual cross- 
sectional area. K is dependent not only upon the optical 
properties of the particle, but also the ratio of particle 

2 Tf 3 

radius (a) to wavelength (x) , in the form a = — — - . There 

are three general regions of scattering determined by the 
value of 

1) Rayleigh scattering region 

2) Mie scattering region 

3) Optical region 

In the Rayleigh region a-*-o, and the particle size is much 
smaller than the wavelength, (an example is scattering from 
air molecules). In this case K increases rapidly with a, 
behaving like at a = o. In the Mie region, the particle 
size becomes comparable with the wavelength and interference 
between the incident light and the scattered light results in 
a value of K which generally oscillates about a mean as a 


2 


increases. In the Optical region the particle is much larger 
than the wavelength and the laws of Geometrical Optics hold. 

The shape of the K vs a scattering curve for forward scattering 
from spherical water droplets is reproduced from reference 4 in 
Figure 1, and is typical in shape of those for most other materials. 


Light which is scattered from one particle may be 
rescattered from a second particle. The cross-sections for 
this multiple scattering are much smaller than those for single 
scattering and will be ignored here. 


For single scattering we can write, using the definition 


of K; 

_dl 

Idx 



a 2 dN' (a) 


o 

where I is the intensity of the light beam, 
x is distance along the beam, and 
dN' (a) is the number of scattering particles of radius 
between a and a + da per unit of length along x. 


( 1 ) 


This equation will serve as the basis for the scattering 
estimates . 


B. Condensation Mechanism 

Liquid and solid condensates, such as water droplets 
or ice crystals can be formed in an unheated air jet when it 
expands to a temperature which is below either the dew point 
of the jet air or the surrounding ambient air. The latter 
case becomes important where the jet air is dried to a very low 
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dew point. Then the ambient, wet air is entrained into the 
jet by turbulent mixing and is cooled by conduction of heat to 
the jet air. At sufficiently low temperatures (high jet 
Mach numbers) the condensed water droplets thus formed may 
freeze into ice crystals. At still lower temperatures, C0 2 
or even 0 2 and N 2 in the air may condense to form solid 
crystals. The dominant scattering specie, therefore is deter- 
mined by the jet operating conditions of stagnation temperature 
and Mach number. 

The dynamics of the formation of solid or liquid con- 
densate particles are extremely complex, since the particles 
do not form instantaneously once the dew or freezing points 
are reached, i.e. some degree of supersaturation always occurs. 
However, if formation is delayed by supersaturation the result 
he fewer particles being formed at a given downstream 
location in the jet. The net result of assuming no supersatura- 
tion then, will always ensure we do not underestimate the total 
number. Another complication occurs in attempting to describe 
the spatial distribution of particles in the turbulent jet. 

The mechanism of mixing by turbulence is a complicated one, 
and in fact its study is one of the basic objectives of crossed- 
beam studies. A further, and again "safe”, simplification will 
therefore be made with respect to mixing of the condensed 
particles. It will be assumed that the particles are uniformly 
mixed throughout that portion of the jet which is below dew 
point. The validity of this assumption improves with increasing 
downstream locations. 
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( 2 ) 


If the beam passes through a length L where the 
temperature is below dew point and the beam diameter is d, 
there is a volume of gas 

V = 

4 

within the beam at any instant of time. If the partial 
pressure of condensible vapor is P (T ) = P before being 

V & Vcl ° 

cooled (temperature T fl ) and is P y (T\) = after being 

cooled, (temperature T, ) then the mole fraction of condensate 
P va P . J T 

lost is p — - where P & is the partial pressure of 

3L &. 1 

air and is assumed equal to the jet (or ambient) pressure. 
Latent heat added by the condensing processes is assumed to 
be small, and not to affect the jet temperature. Also, TL 
should be the average temperature in the jet at the downstream 
location of the beam. The ratio of total mass of condensible 
vapor to total mass of air is therefore 


m 


v 


M 


v 


m_ 


M. 


v_a 

a 


1 - 


P . T 
vj L a 

P T7 
va j 


( 3 ) 


where M y , M a are the molecular weight of vapor and air 
P v j = P g j is the saturation partial pressure at 
temperature T. . 

Having calculated the total mass of condensible vapor 
it is necessary, in order to determine scattering effects to 
estimate particle size distribution. The common probability 
distribution functions (e.g. Gaussian) will not suffice here 
since they allow a finite number to exist at zero and negative 
sizes. A well behaved distribution for our purposes is a 
Gamma distribution and this is discussed in the following 
section. 
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C. Particle Size Distribution 


The Gamma probability distribution has the form^^ 


/ \ a n “itia 
p(a) - A a e 


( 4 ) 


There are several conditions placed upon equation 
(4) which determine the parameters A, n, m. These are 


/ 


1. I p(a) da = 1 
o 


2 . 


00 

/ 


p(a) a da - a^ = mean particle radius 


3. 


J P< a > < a - a m ) 2 da =<(a - a m ) 2 >= 


<72 

2 . 


, i < (a - a ) S 

The shape variable a' is defined by — = m 


cf 2 


m 


Using these conditions we obtain 

( CT/af 

1 . A = — m 


r(c ) 


2. m = 


<T 


m 


3. 


n = 


£T - 1 


or 



1 

N 


d N(a) 

Hi" 


(5) 


Figure 2 show some plots of p(a) for <T = 4, 8, 16, 24, 36, 
oo and a^ - 2(i. In the limit cy*- <» equation (5) defines a 
3 - functional, the particular case where all particles are of 

the same size a m . Since p(a) da = ^ we see that the mean 


6 



particle size and dispersion of particle sizes are easily 
controlled independently through a m and a r . Using equations 
2, 3 and 5 we can now estimate the total number of particles 
and their distribution. 


The volume of each particle of condensate is 4 rra^ 
.4 3 

and its mass is ira p c> The total mass m c of condensate of 
N particles distributed about a mean size according to 


equation 5 then is 

, CO 

r 


00 


m c / 3 ^ P c a 3 dN = 7 t p c N j a^ p(a) da 


oo 


3 ir N p c fTTFT 


<r -° 1 


a 


a 


m 


a 3 a 
a 

m 


<T a 
a 

m 


da 


4 T a 3 p LS. ± m 3L+ 
3 m 1 C ZT2 


( 6 ) 


Now m c - m v , the mass of vapor which could be present without 
condensation, and letting m a = P a V = p a rrd 2 L we obtain from 
equations (3) and (6) , 


N 



M v 




2 

<r 

(cr + 1) { cH- 2) 

(7) 


Here we have assumed the jet saturated so that p . — p 

vj sa 

and that the relative humidity of the ambient air va ^ 

P sa 

Equations (5) and (7) thus allow the determination of particle 
numbers and sizes. It now remains only to calculate the 
scattered signal using equation (1). 
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u. Scattering Signal Determination 


Using Equation (1) we can integrate to obtain 


in I = 


00 

-/ 


f 2 

j K(a) Tra^ d N' (a) dx 


o o 


If a is taken independent of x, i.e. particles are 
uniformly distributed within the beam, then 

in I - - v [k( cl) a 2 d N' (a) l 


-I' 


- if | K(a) a d N(a) 

r 


= - 7 r 


N J 


K(a) a p(a) da 


or 


- tt K eff (a m ) a m N 


nr m 


I_ 

I. 


= ex P - ^ eff (a m ) a m N 


00 

/ 


2 .-I 

J 


(a — I ^ 

\ 3 m / 


( 8 ) 

( 9 ) 


where K gff (c^) = j K(a) 

o 

We note that for <r - p(a) = £ (a - a ) 

V 

and K eff = K(a m ) . 

The total number of particles is a function of the 

dispersion in particle sizes through the factor ( + 1) ( <T + 2) 

6 2 
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F° r CT - <x) j this factor is unity, increasing to 6 for (f - 1. 


In review then, the pertinent equations for the deter- 
mination of light extinction by scattering from condensates 
are as follows: 



-f K eff 


Q-) a 


m 


) 

N i 


( 8 ) 







P 


P . T \ 
sj a 

P~^ t: 

sa j / 



(T 2 

( CT + 1) ( (T + 2) 

( 10 ) 



Restricted by the following assumptions: 

1. No supersaturation 

2. Uniform spatial distribution of particles 

3. Gamma distribution of particle sizes 

4. Uniform temperature in jet where saturation occurs. 
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Ill* COMPUTATION OF SCATTERING SIGNAL 


A. Computer Program 

Equations 5, 8, 9 and 10 comprise the equations necessary 
to calculate expected scattering through a jet with condensing 
vapors. They apply obviously only to mean scattering, considera- 
tion of the turbulent statistics is left to a following section 
The solution to the equations is straight-forward using numerical 
techniques, and a program has been written in Fortran IV for the 
determination of scattering levels. 

Input constants required are as follows: 

d beam diameter 

L effective path length within the jet 

X wavelength 

^ mean particle size 

cf inverse of dispersion in particle size 

molecular weight of vapor 
M a molecular weight of air 

p c j density of condensate at jet temperature 

p a density of condensate at air temperature 

(3 relative humidity of inlet air 

P a pressure of ambient air 

P ga vapor saturation pressure at ambient 
air temperature 

p g j vapor saturation pressure at jet 
temperature 

K(a) scattering cross-section 
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The material properties are all obtainable from handbooks 
(e.g. reference 6). The scattering cross-section K(a) can 
be found tabulated in a variety of texts (e.g. reference 7). 
The humidity of the inlet air (3 must be measured or at least 
estimated. This leaves the particle sizes, as characterized 
by and & as the main unknown inputs. Only educated 
guesses can be made of the likely ranges of these parameters 
and then, by performing calculations, the sensitivity of the 
results to variations in the parameters can be determined. 
Such calculations are reported in a later section. 
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B. Limiting Solutions 


Although in general, the integral for K gff in equation 9 
cannot be solved directly and we roust resort to numerical 
solution, there are particular cases where direct, closed- form 
solutions can be found. Of these, two limiting cases are of 
importance to us, viz. , the Rayleigh region and the Optical 
region, where we can substitute analytical froms for K(a). 

1 • Optical Region 


ine simplest of the 




K ( a) = K = const, 
o 


Then equation 9 becomes 


K 


eff 


CO 

■ / ft) 


p(a) da 


Using equation 5, for p(a) this becomes 

cr+ l 


K eff = K o 


<r 


(ii) 


Rather than give the complete form for AI here, we will 

determine the ratio for a general value of ry to that for 
— U ‘ ^ 1 < 3 ~ 

CT - that is . Further let us assume (although this 

is is not necessary, and is done only for convenience) that the 
exponent of equation 8 is sufficiently small that we can write 

2 „ 

( 12 ) 


AI v 

r~ = t K a 
1 Q eff m 


N 


Then 


AI 

AT 


< K effV 

< K eff)6 


N 


N, 


(13) 
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From equation 11, 

^ K effV 

and from equation 10 


<T+ 1 


N C 

W, 


( <? + 2) ( cr+ l) 


( 14 ) 


(15) 


Using equations 14 and 15 with 13 


A I, 


aT 


cr 

cs* + 2 


(16) 


/j X ry 

As an example, for M > cp > 4, l‘> — - — ^ ± 

AI , 3 


representing a rather small error involved in ignoring the 
particle size distribution. 


2. 


Rayleigh Region 

In the Rayleigh region, as a 


o, K(a) r^a. 


Thus we may write K(a,A) = K(a m ) — ^ 


m 


Then equation 9 becomes 


K - K(a ) 
eff m 


CO 

/( 


a 


m 


p(a) da 


o 


Again, using equation 5, this becomes 


K eff ~ K(a m ) 


( CT + 5) ( <T ± 4)( <r + 3) ( CT + 2) ( g- + 1 ) 

a 73 


( 17 ) 
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Proceeding as before, 


AV 

^7 


( <r+ 5)( o~ + 4) ( (T4- 3") 
<T 3 


In this case, for oo )> (T )> 4,1 


AV 

"aTT 


< 


7.87 


( 18 ) 


We have a sizeable possible error for this case and a 
complete reversal of it from the previous case; i.e. ignoring 
the distribution of particle sizes underestimates rather than 
overestimates the scattering, in contradistinction to the 
observed behaior in the optical region. 


In general, K(oc) in the Hie region, is very complicated 
in shape, and a particular value of a m and CP may result in 
particle sizes lying within two regions. Thus, although these 
limiting cases are of interest and indicate trends with a 

m 

and C T we need numerical solutions to study the complete behavior. 


C . A Sample Problem 

Of particular interest to crossed beam measurements 
is the scattering to be expected from water vapor condensing 
within a turbulent jet at a wavelength near 4.3 p. The jet 
velocity is obtained by expanding from a compressed air supply, 
which in turn has a very low absolute humidity since most of 
the water vapor condenses out during interstage cooling. 

Water vapor condensation therefore occurs when ambient air is 
mixed into the cool jet, and we need to know the ambient 
humidity. For an ambient temperature of 27°C realistic room 
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dew point is 15°C which corresponds to p s .5 at 760 torr 

( 8 ) 

pressure. Crossed-beam studies at IIT Research Institute 
are made using a 2.54 cm jet at a typical flow Mach number of 
M = .62. The jet temperature, expanding from an ambient 
temperature of 27°C is then 4°C. Thus water droplets, not 
ice crystals are formed. The following conditions were chosen 
for sample calculations: 


d 

L 

X 

M 

v 

M„ 


fJ cj 

Pa 

V 

P a 

P sa 

P sj 

K(a) 

cr 


m 


. 1cm 
10cm 

.2 to 20p, 

18 

29 

1.000 gm/cc at Tj = 4°C 
.001177 gm/cc at T - 27°C 
0.5 

760 torr 

26.739 torr at T = 27°C 

cl 

6.101 torr at Tj = 4°C 
taken from Figure 1. 

4 , 16 , 24 , oo 
1 to 20ji 


The results for = <», a m = lb- are shown in Figure 3 

plotted as vs <x . Since < 107,, the exponential of 

i o i o 

equation 8 may be simplified by the assumption e"^ = 1 - y and 

^ — from equation 8 and 10. Thus the results of Figure 3 
1 o a m 

may be readily extrapolated to other values of a . 
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The effect of shape of the distribution function (T ) 
can be seen from Figure 4. The curves have been normalized by 
the value at S' - oo . Although ? has a pronounced effect at 
small a as noted earlier, the absolute value of the scattering 
drops off rapidly below a= 1 (c.f. Figure 3). The result is 
that we are concerned primarily with a >l and here the ratio 

between 3 and 2/3 for S' ^>4, As a first approximation 

6 

then, we can ignore the correction for <r and after determining 
the conditions for maximum scattering, reassess the effect 
°f ** * I n Figure 5,-^- is plotted for two values of X : one 
corresponds to the CO2 band center at 4.3 \± and the other corres- 
ponds to the ultra violet oxygen bands (Schumann- Range) at ,18[i. 
The figure points up clearly the increased effect of peak scatter- 
ing at small wavelengths. The peak scattering occurs at a » # 15u 
or a = 6 for X - .1% and at = 3.4p. or a - 5 for X = 4.3p. 

At these values of qt , the correction for from Figure 3 is 
less than 40/ o for 4<^'T<' 00 . At X = 4.3p the maximum scattering 
expected for our test case then is about 3% and at x - . 18p is 

65%. Crossed-beam studies are made with ^ 50% to give 

Io 

an optimum signal- to-noise ratio. Hence at x = 4.3p scattering 
from water droplets should be negligible but not necessarily so 
at X — . 18p„ . This points up the advantage of performing the 
crossed-beam studies using the infra-red system. 

In order to apply these results to experiments at condi- 
tions other than those assumed here, a set of extrapolation rules 
has been prepared and these are discussed in the following 
section. 
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D. Extrapolation to Other Flow Conditions 

As observed previously, yi < 10% snd we can approxi- 

l o 

mate the exponential by the assumption e ^ = 1 - y. Then 

equations 8 and 10 indicate a linear relation between yi and 

, L, and H = (p - . Extrapolation to 

a m Pcj P a ^sa 1 ^ 

conditions other than those of the sample calculation then 
becomes straightforward. 

a) To correct for particle size the results of Figure 3 
need only be divided by the mean particle radius in microns. 

b) To correct for path length L, the results of Figure 3 
should be multiplied by L/10 with L in cm. 

c) The density of condensate is p . = 1 gm/cc to within 2% 

C J 

accuracy from 0°C to 70°C . Consequently the main correction in 
- 2 - is for p . Since p a = .001177 gm/cc (at T a = 300°K and 

Pcj a a 

p = 760 torr) was used in the sample calculation we can correct 
for ambient air density by multiplying the results of Figure 3 

by = .395 ^ , where P a is in torr, T & in °K. 


d) The remaining correc tion ' is really a correction for 

T a and Tj . Thus to correct for humidity effects, H can be 
calculated from P and P • (determined from Reference 6 and 

S3 S J 

plotted as in Figure 6) and the known value of p. Since 

_ O 

H = 9.574 x 10 for the sample calculation, the results of 

IJ 

1 

Figure 3 should then be multiplied by — . 

6 9.564 x 10-3 
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IV. EXPERIMENTAL STUDIES 


Single beam autocovariance measurements were taken 
using the beam-splitter arrangement on the infra-red system, 
at a downstream location of 4 jet diameters. The beam-splitter 
arrangement gives a measurement independent of electronic and 
detector noise. The reading obtained is of course the mean- 
square fluctuating signal and we are interested in comparisons 
with the extinction of the mean signal. To accomplish such a 
comparison we must make some assumptions concerning the statis- 
tical behavior and interdependency of the scattering and ab- 
sorbing signals. 


A * Relation of Measured Signals to Mean Extinction 

It has already been assumed that the scattering particles 
are completely mixed within the jet by the action of turbulence. 
If it is assumed that density gradients are mixed in the same 
way, we would expect the statistical behavior of scattering and 

absorption to be identical and we will make such an assumption 
here . 


Let V *a be the instantaneous a.c. signals from 
scattering and absorption and AI g and Al a be the corresponding 
mean signal extinctions 


then i = i + i 
s a 


(19) 


and 


<i 2 >= <<i s + i a ) 2 >- <is> + < 1 a > + < i s V < 20 > 
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Since the mixing processes are taken identical for scattering 
and absorption 


Then 


/h 

S 

-a ? > 

< i s 1 .> 


Ala 

AI 8 

< i s> . 



<i 2 > 


+ AIq + AI a AI S 


from which 



Al a 



2 


s 


( 21 ) 


( 22 ) 


(23) 


The + sign holds where scattering dominates (in wings of band) 
whereas the - sign holds where absorption dominates (in band 
center). 

2 

The ratio s is measurable with the infra-red system 
by making measurements in the 4.3 p band and obtaining <i 2 > 
and then taking measurements off-band where only scattering 
should exist, and obtaining <( ig)>» Early results with the 
infra-red system showed evidence of oil particles causing 
scattering but this problem has since been alleviated and it 
is believed, is no longer a problem. The total absorption 
extinction Al a is measured by taking d.c. intensity readings 
across the 4.3p band and correcting for that portion of the 
path which does not pass through the jet. It is assumed that 
Al g « Al a and can be neglected in this measurement. This 
limits us to the - sign solution of equation (23). Such measure- 
ments are reported in the next section. 
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B, Observations 


2 

Profiles for <( i >,1 and AI q taken across the 4.3^ band 
are reproduced in Figure 7. The estimated value for / i 2 \ 

x S' 

2 

was assumed constant in the calculation of s . The ratio 
AI s 

was calculated from equation 23 (using the minus sign) and 
a AI 

the estimate of — ^ determined (See figure 8). Thus the ex- 

1 C 

periment indicates that the mean extinction due to scattering 

A I 

is about 67o of the light incident intensity i.e. — 2: .06 

x o 

in the center of the band. 


The actual condition for the experiment were 
= 29°C, Tj = 4°C, P = .4, M = .62. The dew point occurs 
at T d = 19°C. If we assume the stagnation temperature constant, 
equal to T g throughout the jet, then we can use the local 
isentropic relation to calculate the local speed. 


273 + T 

i.e. 2 

273 + T 

d 


302 

292 


1+0.2M 2 


2 

I U(y/D)\ 

U(o) / 


where U(y/D) is the speed at the radial position y/D where 

T = T d and U(o) is the center- line speed. This gives ^(y/D) _ 

, /ii\ U(o) 

.62, which occurs v > at y - D/2. Thus the width of 

jet of interest, L = 2y = D = 2.54 cm, instead of L = 10 cm as 
used in Figure 3 and we must correct the results (using extrapol- 
ation procedure b) of Section III D) by a factor of ~ 
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The correction for ZSl is less than 1 % and will be 

ignored. However a correction is required for H in order to 

use Figure 3. From Figure 6 } taking T = 29°C ) T. = 19°C, 

P P . 3 3 

we obtain „- S - a = .04 and = .021 then H = ,04(.4 - .02' 


.021 then H = ,04(.4 - .021 x 


309 a _3 a 

= «0076 = 7.6 x 10 . The correction (using the extrapol- 
ation precedure c) of Section III D) then is - — X ^ - ■— 0; 0.8 

08 9.6 x 10-3 

and the total correction is 7 - - = .2. The experimental 


measurement indicates that 


A ^ Q 

=— is about .06, under the 
*o 


condition tested. If the conditions had been those of the sample 

A! s 06 

problem, we would have obtained a measured — = ’ — = .30. 

o . 2 

This number is much higher than the maximum value s = 

shown in Figure 5 at 4.3p and would lead us to the I o 


conclusion that water vapour condensate is not sufficiently 
large to account for the off-band signal. 
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V. DISCUSSION 


The analysis presented here, designed to give a maximum 
estimate., predicts that scattering from condensed water droplets 
should be no problem at X = 4.3p,. Experimentally, however, we 
find that scattering levels, although not dominant, are an 
appreciable fraction of the total signal. Although the estimates, 
both analytical and experimental are rough, accuracy is suffi- 
ciently good that we must conclude that there is scattering 
occuring from particles other than water droplets. Earlier 
studies showed that oil droplets in the IITRI facility domin- 
ated the scattering signal, and although this has been reduced, 
it would appear that it still amounts to at least an order of 
magnitude more signal than can be expected from water vapor 
c ondensate . 
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Forward Scattering Cross-Section, 



. J L i ...... I i _1 J 

.2 .4 .6 .8 1 2 4 



Figure 1 CROSS SECTIONS FOR FORWARD SCATTERING 
FROM SPHERICAL WATER DROPLETS. 
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Probability Density, p (a) (micron) 


l I 



Droplet Radius, a microns 


Figure 2 GAMMA PROBABILITY DISTRIBUTION FOR MEAN 

DROPLET RADIUS a - 1 MICRON. 

m 
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Mean Square Signal 

Mean Signals, Arbitrary Units Arbitrary Units 




Figure 7 MEASURED MEAN AND FLUCTUATING SIGNALS 
FROM I.R. SYSTEM. 


29 



Scattering Extinction x 100% 



Wavelength, microns 
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